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Abstract 

The Green-Kubo relation, the Einstein relation, and the fluctuation-response re- 
lation are representative universal relations among measurable quantities that are 
valid in the linear response regime. We provide pedagogical proofs of these universal 
relations for stochastic many-body systems. Through these simple proofs, we char- 
acterize the three relations as follows. The Green-Kubo relation is a direct result 
of the local detailed balance condition, the fluctuation-response relation represents 
the dynamic extension of both the Green-Kubo relation and the fluctuation relation 
in equilibrium statistical mechanics, and the Einstein relation can be understood 
by considering thermodynamics. We also clarify the interrelationships among the 
universal relations. 
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1 Introduction 



Let us consider the general class of systems exhibiting the macroscopic flow 
of some quantity. Familiar examples include electric conduction systems, heat 
conduction systems and driven colloidal systems. Such systems often exist in 
states in which the current is stationary. We refer to these states as nonequi- 
lihrium steady states. Nonequilibrium steady states are the simplest type of 
nonequilibrium states. However, even for the class consisting of such simple 
states, there presently exists no universal framework to calculate macroscopic 
properties of systems from microscopic information, in contrast to the class of 
equilibrium states, for which macroscopic properties can be calculated from 
the partition function of the system, employing the well-established theory of 
statistical mechanics. 
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At the present time, there exists a theory of nonequihbrium systems only for 
the regime very near equihbrium. This so-called linear response theory con- 
sists of a set of universal relations, which includes the Green-Kubo relation, the 
fluctuation-response relation and the Einstein-relation. These relations relate 
the dynamical properties of fluctuations in systems under equilibrium condi- 
tions to linear transport properties of nonequihbrium systems. The validity 
of these relations has been confirmed by both experimental and theoretical 
investigations [1]. 

In recent years, properties regarding fluctuations and linear responses to per- 
turbations have been investigated for nonequihbrium states even outside the 
linear response regime, in steady state systems [2], aging systems [3,4,5] and 
other systems [6,7]. Obviously, we cannot expect the relations proved within 
linear response theory to be generally valid outside the linear response regime. 
However, there have been several relations proposed and investigated recently 
that represent extensions of the relations in linear response theory to systems 
far from equilibrium [8,9,10,11,12,13,14]. 

We wish to acquire a systematic understanding of these extended relations 
among fluctuations and response properties, aiming at the construction of a 
statistical mechanics of nonequihbrium systems far from equilibrium. For this 
purpose, first it is necessary to grasp the essence of the basic framework of 
linear response theory, in particular focusing on the elucidation of the interre- 
lations among the linear response relations, thermodynamics and equilibrium 
statistical mechanics. 

Of course, all of the relations in linear response theory can be derived by 
applying a formal pcrturbativc expansion to quantum or classical mechanical 
descriptions on the basis of equilibrium statistical mechanics [15,16,17,18,19]. 
However, such formal derivations of these relations does not illuminate the 
physical picture (see Section 8). In this paper, instead of applying a formal 
perturbation theory to quantum or classical mechanical systems, we attempt 
to elucidate the essence of linear response theory by studying stochastic many- 
body systems. 

Among the various types of stochastic many-body systems, we choose to in- 
vestigate two kinds of one-dimensional nonequihbrium lattice gases in which 
particles hop stochastically from one lattice site to another according to a 
given transition rate [20,21,22,23]. Specifically, we study a driven lattice gas, 
in which particles are driven by an external driving force in the bulk, and a 
boundary-driven lattice gas, in which particles are driven by the difference in 
chemical potential between two particle reservoirs in contact with the bound- 
aries. Although we choose these models for their mathematical simplicity, it is 
straightforward to apply the method of proof to other stochastic many-body 
systems, such as Langevin systems that describe the motion of many beads in 
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a liquid. 



This paper is organized as follows. In Section 2, we state our physical point 
of view and the basic assumptions used in the derivation of the linear re- 
sponse relations. In Section 3, we introduce a one-dimensional driven lattice 
gas. In Section 4, wc present a simple derivation of the Green-Kubo relation, 
employing a key identity, with which various universal relations characteriz- 
ing nonequilibrium steady states can be derived. In Section 5, we derive the 
fluctuation-response relations and elucidate its connection to equilibrium sta- 
tistical mechanics. In Section 6, we derive the Einstein relation for the model, 
and clarify the connections among these three universal relations and thermo- 
dynamics. In Section 7, we derive the Green-Kubo relation for a boundary- 
driven lattice gas, which is a different type of nonequilibrium lattice gas, in 
order to demonstrate the generality of our argument. Finally, in Section 8, we 
give several important remarks concerning linear response theory. 



2 Basic assumption 

Let us recall that the macroscopic properties of an equilibrium system can be 
investigated with statistical mechanics using an effective multi-dimensional 
variable o; and the corresponding effective Hamiltonian H{(x). More precisely, 
in order to obtain such a description, the variable a, described statistical 
mechanically, need not constitute a complete set of mechanical variables that 
specifies the microscopic state in classical or quantum systems. For example, 
consider a system consisting of many beads suspended in a liquid. In this case, 
oc can be regarded as representing the set of centers of mass of the beads, 
with H{a) describing the effective interactions between them. The effect of 
the liquid can be included in this effective Hamiltonian H{cx). 

Here, we suppose that the set of variables (ck, F) provides a complete mechan- 
ical description of the system under consideration. For the bead system, F is 
defined so that (a, F) represents all the positions and conjugate momenta of 
the atoms composing the beads and the liquid in an adiabatic container. 

Now, we consider a nonequilibrium system which is such that the degrees of 
freedom represented by F remain in equilibrium, while only those represented 
by cx are in a nonequilibrium state. It is not certain that such an assumption 
would be valid, and in general the degree of error it introduces would depend 
on the system we consider. Physically, this assumption seems reasonable in 
the case that the typical time scale of a. Ta, is much longer than that of 
F, rr. We consider such and introduce a time interval At satisfying 

Tr -C At -C Tq. In the argument below, we refer the subsystem composed of 
the degrees of freedom F as the heat bath. Because this heat bath is assumed 
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to be always in equilibrium, its temperature, T, is defined for nonequilibrium 
states of the remaining degree of freedom. 



In order to describe the nonequilibrium states of such systems, we need to 
determine the time evolution of the variable ot on the time scale of the interval 

At. Due to the separation of time scales represented by the relation rr ^ At, 
it is reasonable to assume that a is described by a Markov stochastic process 
with a transition probability T(a cx'). 

The fundamental assumption in this paper is expressed by the local detailed 
balance condition 

T{a ^ oc') _ 0Qi^^^') , . 

T{cx> ^ a) " 

where (3 — 1/T. Here, the Boltzmann constant is set to unity. The quantity 
Q{cx ol') represents the energy transferred from this heat bath, which we 
call heat, during the time evolution symbolized by a — > a' taking place over 
the time interval At. 

In the system under equilibrium conditions, because the energy conservation 
expressed as Q{cx ol') = H{ol') — H{a) should hold, (1) becomes the 
detailed balance condition of the transition probability: 

p"^°(a)r(a ^ a') = p'^''{a')T{a' a), (2) 

with respect to the canonical distribution 

p-f3H{c.) 

r'\cx) = (3) 

where Z is a normalization constant. Combining (2) with the normalization 
condition of the probabihty, J2a' T{cx. — > cx!) = 1, we derive 

p^'^'^(a) = ^p^'^"(a')T(a' ^ ot). (4) 

a' 

This equality implies that (3) is the stationary distribution of ol. This result 
is consistent with equilibrium statistical mechanics. Also, it turns out that 
the detailed balance condition (2) implies time-reversal symmetry in the sense 
that the two-time joint probability p'^^'^(a)T(a — > ol') is invariant under the 
exchange of ot and a'. 

By contrast, in the case of a nonequilibrium state, it is not easy to obtain a 
physical interpretation of (1). We explain the origin of the name local detailed 
balance condition in Section 3, and demonstrate that all the universal relations 
can be derived from this assumption in Sections 4, 5, and 6. In the remainder 
of this section, we explain how the local detailed balance condition (1) places 
restrictions on the allowed types of interactions with the heat bath. 
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In order to demonstrate this point explicitly, we consider a many-body bead 
system in which an external driving force is applied to the beads. We wish 
to derive the transition probability T{{rj} — > {r'j}) over an infinitesimally 
small time interval At, using the local detailed balance condition, where 
(i = 1, • • • ,N) represents the position of the i-th bead. 

First, note that the mechanical force acting on the i-th bead is given by 

where Uq is a one-body potential (e.g. a spatially periodic potential), U is 
the interaction potential between beads, and / is an external driving force. 
Then, because the change in momentum is much smaller than the dissipative 
force in this situation, the averaged motion caused by the force Fi[{rj}) is 
determined by the force balance equation 

-7((r-:)-(O) + F,Ai = 0, (6) 

where 7 represents the dissipation constant, and Fi is the average of Fi during 
the transition, given by 

Fi = l[Fi{{r,{t)}) + F,{{r,{t + At)})]. (7) 



Under the influence of the interactions with molecules in the liquid, the posi- 
tions of the beads fluctuate around this averaged motion. Because we choose 
At to be much larger than the time scale characterizing the dynamics of the 
molecules, we expect that these fluctuations exhibit a Gaussian distribution 
with a dispersion proportional to At~^. That is, it is reasonable to express the 
fluctuations by the transition probability 

Tiivj} - {r;.}) = le-^*iHELl7^.-F.P^ (8) 

where y is a normalization factor, 6 is a parameter which represents the 
strength of the interaction with the heat bath, and Vi should be interpreted 
as 

We wish to determine b in (8) from the local detailed balance condition. In 
order to do this, we consider the following ratio of transition probabilities: 

In this expression, it should be noted that the quantity At J2iLi i^iFi is the 
work done by the force Fi, which is equal to the heat dissipated into the heat 
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bath [24] during the transition from {vi} to {r-}. Then, by comparing (10) 
with the local detailed balance condition (1), we find the relation 

b = 7T. (11) 

In this way, we have determined the transition probability (8), with the explicit 
form of b from the local detailed balance condition (1). 

Note that the stochastic process given by the transition probabihty (8) can be 
expressed as 

j{r[ r,) = FAt + Aiy„ (12) 
where APVj represents Brownian motion satisfying 

{l^WAWj) = 2b5ijM. (13) 

Then, taking the limit At — > 0, we obtain the Langevin equation 

where represents zero-mean Gaussian white noise with 

{mut'))-'^b^i3m-t'). (15) 

Here, / is a unit matrix. The condition (11) is called the fluctuation- dissipation 
relation of the second kind for the Langevin equation (14) with (15). 



3 Model 



We can derive all the universal relations in linear response theory for the 
Langevin equation (14) with (11) and (15). However, this derivation is compli- 
cated when mathematical rigor is maintained. Instead of the Langevin model, 
we study a lattice gas, which also serves as a model of many-body beads sys- 
tems, to avoid the complicated analysis that is irrelevant to the essence of 
linear response theory. 

In a lattice gas, particles hop stochastically from one lattice site to another 
at a given transition rate. Particularly, when two or more particles cannot be 
placed at one site, the equilibrium state is described by a Hamiltonian which 
is equivalent to that for a spin model with a variable transformation. For 
mathematical simplicity, the form of the Hamiltonian usually corresponds to 
the Ising model. In Subsection 3.1, we explain the one-dimensional equilibrium 
lattice gas. A two or higher dimensional model can be easily interpreted from 
the one-dimensional case. 
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A driven lattice gas, introduced in Subsection 3.2, is defined by minimally 
modifying to the equilibrium lattice gas so that it exhibits nonequilibrium 
steady states. The model was numerically and theoretically investigated in 
Refs. [21,22,23]. 

3.1 One- dimensional equilibrium lattice gas 

We first define an occupation variable, r]^, on each site a; in a one-dimensional 
periodic lattice x = 1,2, ■ ■ ■ , L. r]^ = I when the site x is occupied by a particle 
and Tjx = if X is unoccupied. A periodic boundary condition is imposed by 
introducing a boundary site ai x — and setting r)o — rji. The collection of 
all occupation variables {fji,--- ,r]i) represents the particle positions in the 
system. We denote [rji, - ■ ■ , rji) as rj, which is referred to as the configuration. 
For convenience, we use r]^ to represent the configuration when the value of rj^ 
is replaced with that of l—f]x- In a similar way, rj^^ represents the configuration 
when the values of rj^ and r)y are replaced by 1 — t^^; and l—fiy. (Note that the 
configuration after exchanging the value of rj^ and rjy is represented by this 
T]^^.) The interaction between the particles is described by the Hamiltonian 

L-l 

Hoiri) = - 1^ VxVx+i- (16) 

x=0 

The time evolution of rj is expressed as follows: At each time step, a nearest 
neighbor pair x and y, < x,y < L, is randomly selected and the values of rjx 
and rjy are exchanged using the exchange probability Co{x, y; rj) that satisfies 

co{x, y; m) = co(x, y; r/^^) e-^«°(^^''^^), (17) 

with 

Qoiri^r)') = Hoir}')-Ho{ri). (18) 
Qoiv ~^ V^^) corresponds to the heat absorbed from the heat bath for the 
configuration change 77 —>■ rj^'^. We also define CQ{x,y; 77) = for pairs x and 
y such that \x — y\ ^ 1. We regard this step as a Monte Carlo step. 

Examples of common concrete forms of cq{x, y\ r\) are the Metropolis method 

Co(x, y; 77) = min{l, g-'^^o^^^^'')}, (19) 
the heat bath method 

and the exponential rule 

Co(x, y; T7) = const. Q-^Q^kn^^^) ^ (21) 
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where const, is chosen so as to satisfy the condition that co(x,y;ri) < 1. It 
can be easily checked that (19), (20) and (21) satisfy (17). 

The transition probabihty To(77 — > t]') for one Monte Carlo step is written as 

To(r]^r]')^^co(x,y;rj) (22) 

for Tj' = rj^y 7^ 77, and To{t] r/) is determined by the normahzation condition 
of the probability 

^Toirj ^ ri') = 1. (23) 
v' 

Note that in the argument, To(r7 rj') is different from co{x,y; rj); the tran- 
sition probability, To(ri — >■ rj'), is a function of two configurations, while the 
exchange probability, co{x, y; 77), is a function of the sites x, y and the config- 
uration 77. 

From (17) and (22), we find that To(77 77') satisfies the local detailed balance 
condition 

Toirj V') ^ pQo^rj^rj') (2d) 

Here, using (18) in this equilibrium case, the condition (24) becomes the de- 
tailed balance condition 

p^^^{ri)Uri ^ 77') = ^^^-^(770^0(77' ^ 77), (25) 

with respect to the canonical distribution 

-/3ifo(r?) L 

p'^^rj) = ^^5(E V.-N). (26) 
Here, Z is a normalization constant, and is the total number of particles, 

N^j^Vx, (27) 

x=l 

which is conserved through this time evolution, p = N/ L is introduced to be a 
parameter of the model. Similar to Section 2, using (25), we can confirm that 
(26) is the stationary distribution of this model in the equilibrium case. 



3.2 One- dimensional driven lattice gas 



In this subsection, a driven lattice gas that describes a physical situation where 
the particles are driven by a uniform driving force E in the bulk is introduced 
as a simple extended model of the equilibrium lattice gas. First, we define the 
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net number of particles that hop from x to x + 1 in the configuration change 
r/ — > r/' as 

^ V) = rj,{l - v'M+iCi- - Vx+i) - r/i(l - V.)v.+i{'^ - vUi)- (28) 

Note that ^xiv — ^ "H') takes a value in {±1,0}. Then, the heat Qq given in 
(18) is replaced with 

Q{rj ^ 17') = Hoir^') - Ho{r)) - E ^ *-(^ ^ V), (29) 

x=0 

to consider the effect of the driving force on the time evolution rule, because 
Q{r) — >■ T]') is interpreted as the heat absorbed from the heat bath for the 
configuration change ry — > 77' in one Monte Carlo step. Using this replacement, 
Cq{x, y; rf) in (17) and To(r7 rj') in (22) are replaced by c{x, y; rj) and T{r} 
T]'), respectively. Then, T{r} — * rj') satisfies the local detailed balance condition 

Here, one may consider whether Q{r] ^ t]') given in (29) can be expressed as 
Q{ri ^ V) = E h,{ri') - j:hx{v), (31) 

x=0 x=0 

with a local Hamiltonian 

hxiv) = -VxVx+i - Exf]^. (32) 

If (31) and (32) were valid, this should correspond to the detailed balance 
condition, because (30) can be written as 

T(t7 ^ 77') e-^S^=o = T(t7' ^ 77) e-^Sx=o . (33) 

However, (31) is valid only when rj' = 77^^+^ with < x < L — 2, and is 
invalid when x = L — 1, unless E = 0. More generally, when E ^ 0, (31) 
can be satisfied for the configuration changes with one exception, but not for 
all exchanges, even if we introduce another local Hamiltonian. Due to this 
property, the condition (30) is called the local detailed balance condition. As 
shown in the above argument, the local detailed balance condition does not 
necessarily imply the detailed balance condition, except for the equilibrium 
case. 

Throughout this paper, we regard L Monte Carlo steps as the unit of time, 
which is the number of Monte Carlo steps necessary for all sites to be se- 
lected once on the average. This unit of time is called one MCS (Monte Carlo 
step per site). Then, we denote the history of the configurations from to r 
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MCS as [77] = (77(0), r]{l), ■ ■ ■ , ri^Lr)). The probability of a history [77] in the 
nonequilibrium steady state, -Pg ([77]), is written as 

PUH) = P^ivmnm - Vil)) ■ ■■nniLr - l) ^ rjiLr)), (34) 

where p%{ri) is the stationary distribution in the system with the driving force 
E. Specifically, when £^ = 0, Pe=o,a?"(^) identical to ^'^^'^(77) defined in (26). 
Then, the statistical average of a history dependent quantity A in the steady 
state is denoted as 

{A)%^Y.PE{[ri\)A{[ri\). (35) 



4 Green-Kubo relation 



In this section, we provide a simple proof of the Green-Kubo relation for the 
one-dimensional driven lattice gas. During the derivation, we obtain the key 
identity given in (44). Using this key identity, we can derive the nonlinear 
response relation, the fluctuation theorem and the steady state distribution 
as well as the Green-Kubo relation. 



Main claim 



We first define the spatially averaged current per one MCS at time t, j{t). 
Because the number of particles moving from the site x to the site x-\-l from 
i to i -I- 1 is expressed as Z^^Zg ^xivi^t + k) ^ vi^t + k + j{t) is written 
as 

j{t) = y E E ^MLt + k)^ ri{Lt + k + l)). (36) 

^ x=Cl k=0 

Using the steady state current, J — the conductivity, a, is defined as 

Furthermore, for the r-averaged current 

Jr{[ri])^lj:m, (38) 
its fluctuation is characterized by 
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Then, the Green-Kubo relation connects the conductivity, o", with the current 
fluctuation, by 

^=1- (40) 



Proof: We prove (40) using a rather special situation. Suppose that the statis- 
tical distribution of the configuration 77 at t = is canonical, p'^^^{tj), and the 
driving force E is turned on at t = 0. In this transient process, the probability 
of a history [r]] is given by 

^I(W) ^p'-'^ivmnm ^ r7(l)) • --TiniLr - 1) ^ 77(Lr)). (41) 

Then, from (29), (30) and (38), we can derive 

where [rj] represents the time reversed history {r]{LT), r]{LT — 1), • ■ ■ , r7(0)). 

Because the statistical average of a history dependent quantity A by the prob- 
ability -Pg ([^7]) is written as 

{A)I^Y.PeM)A{HI (43) 
M 

the following identity is derived, using (42): 



(/l)^; = Ee/^^-^^^(M)p|([^])^([^]), 

\ri\ 

^/^^-0ErLJ. ^Y^^ (44) 

where ^([^7]) = A ([17]), and we have used JriXv]) = ~'^t([^]) in deriving 
the second line. The expression (44) is the key identity to derive the universal 
relations in linear response theory. We emphasize that this key identity results 
from the local detailed balance condition (30) and that this identity is valid 
even outside the linear response regime. 

Here, Setting A = in (44) , we obtain 

{Jr)%-\{{l-^-'''''^)Jr)'l (45) 

which is called the nonlinear response relation. Yamada and Kawasaki initially 
derived it [25] and then Kawasaki and Gunton generalized it [26]. Finally, by 
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expanding the right hand side of (45) with respect to E, (45) becomes 



J^-E + 0{E'), (46) 

where in the large r limit, the statistical average in this special situation ( 
is equal to the average in the steady state {)%■ This expression provides (40) 
when the hmit ^ is used. 

Remark: One may note that a long time tail of current fluctuation yields an 
anomalous contribution to the Green-Kubo relation. This long time tail is 
mainly due to momentum conservation [27]. In the one- dimensional driven 
lattice gas, the numerical check indicates that the long time tail of the current 
fluctuation only appear when E 0. (See Refs. [28,29] for the case E ^ 0.) 
We conjecture that mode coupling effects are vanishing when E — 0, due to 
the absence of the momentum degrees of freedom. Thus, we cannot discuss a 
compatibility of the long time tail with the Green-Kubo relation within the 
one-dimensional lattice gas. However, it is generally understood that a long 
time tail does not cause the Green-Kubo relation to breakdown when we take 
a long time limit with a fixed system size [30] . 



4-2 Steady state distribution 

We can show that the key identity (44) yields several relations by substituting 
the appropriate quantities into A. Through understanding of the fluctuation 
theorem, originally discovered in Ref. [31], the importance of (44) has been 
recognized [32,33,34,35]. Indeed, by setting ^([^7]) = 6{ETLJr{[i]]) — a), and by 
setting A([r;]) = 1 in (44), we obtain the fluctuation theorem. It is noteworthy 
that the connection of the fluctuation theorem with the nonlinear response 
relation (45), which is also derived from (44), was noted by Crooks [35]. 

Furthermore, as indicated in Ref. [35], setting ^(['fy]) = d{r) — r]{T)) in the key 
identity (44), the steady state distribution function is obtained as 

p%{n) = ^lim (e-^^^^'^ Sin - 77(0)))'; . (47) 

(47) can be rewritten so that the deviation from the equilibrium distribution 
function, p^^"-{ri), is explicit. That is, 

where ( )En{o)=r] I'epresents the average with the initial condition of 77(0) = rj 
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when t — 0. Here, (48) can be also rewritten as 

p|(r7) = p--(r7) }}me^n=^ ;^(-/3^-^)"((^0">^"-(o)=., (49) 

where ( )^'^^(o)=Tj denotes the cummulant. It should be noted that the expres- 
sions of the steady state distribution (48) and (49) are valid even outside the 
linear response regime. 

Particularly, in the linear response regime, (49) is reduced to 

This expression indicates that the deviation of p%{t]) from the equilibrium 
distribution, p^'^'^{rj), is represented by an entropy production. Zubarev and 
McLennan found a similar expression by studying the nonequilibrium steady 
state distribution [36,37]. 

Furthermore, when the contribution of the order E"^ is considered, p%{tj) can 
be expressed as 

fUv) = p-{ri) Jim e-^^^"<^^>--(°)--^^(^^^^)^(("^)^)-oMo)^.+^(^^) . (51) 

It is interesting to see the fact that the existence of the Green-Kubo relation 
guarantees the convergence of the right-hand side of the expression (51) with 
the limit r ^ 00. 



4.3 Reciprocity relation 

To simply demonstrate the reciprocity relation [38,39], let us consider a rather 
artificial model that consists of two types of particles labeled, (X) and (Y). 
(X) particles are driven by Ei while (Y) by E2. In this situation, Q{r] ^ rj') 
defined in (29) becomes 

Qirj - T7') = Ho{ri')-Ho{ri)-E,j:^^{rj ri') - E^J^ ^^Jrj r/'), (52) 

x=0 x=0 

where $^ and $J are the net numbers of particles labeled (X) and (Y) that 
hop from a; to x + 1 for the configuration change rj rj', respectively. Then, 
we define the spatially averaged currents per one MCS at time t, j^{t) and 
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L-1 L-1 

L 

L-1 L-1 



j'^it) ^ T E E ^^(^(^^ + ^) - ^(^^ + k + 1)), (53) 

x=0 k=0 



fii) ^ 7 E E + via + + 1)). (54) 

Using (53) and (54), the r-averaged currents are defined as 



Jrm^-T.fit), (55) 

^ t=0 

j^m^^jifit). (56) 



^ t=0 



When El ^ and £^2 = 0, we derive the key identity in this situation: 



= (e-/^^-^-^? A)l . (57) 



Setting A — Jj and expanding (57) with respect to Ei, we obtain 

^X.^-o = l^^^^^{'rjX,-^^OiE!). (58) 

Here, wc define the conductivity, (T12, which characterizes the transportation 
of type (Y) particles caused by the driving force Ei (only applied to type (X) 
particles), as 

(T12 = lim — (fY^ . (59) 
Then, (58) can be rewritten as 

.:, = i/JrL(j?j7);^„^^^^„. (60, 



Considering the opposite case, when £^2 7^ and Ei — 0, we define the con- 
ductivity, (T21, as 

(721 = lim ifT ■ (61) 
In a similar way as deriving (60), we can also derive 

From the trivial identity 

JrJrT =(jJjrT ' (63) 

^ ^ / Ei=0,E2=0 \ ^ ^ / Ei=0,E2=0 ^ ' 
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we obtain the reciprocity relation between the conductivities as 

<^12 = <^21- 



(64) 



5 Fluctuation-response relation 

The fluctuation-response relation (or the fluctuation-dissipation relation) rep- 
resents a relation that connects a time dependent response to a time dependent 
external force with a time-correlation function in equilibrium. There are two 
types of fluctuation-response relations. One, which is regarded as an exten- 
sion of the Green-Kubo relation, is when the external force is a time-dependent 
driving force. The other is when the external force is generated from a time 
dependent potential and is regarded as an extension of the fluctuation relation 
that can be derived within equilibrium statistical mechanics. 

5.1 Main claim (time- dependent driving force) 

Let us consider the case that the strength of the driving force E{t) changes 
at every MSG. When the strength is sufficiently weak, the time-dependent 
averaged current {j{t)) (with an initial condition given ait = — oo) is expressed 

by 

oo 

{m) = Y.mE{t-s) + 0{E'), (65) 

s=0 

where R{s) is independent of E, and is called a response function. Note that, 
from the causality, R{t) = for t < 0. In a special situation when the driving 
force E is turned on at t = 0; that is, E{t) — £^(=const.) for t > and 
E{t) = for t < 0, the relation (65) becomes 

{m)% = EY:R{s) + 0{E'). (66) 

On the other hand, setting A = j(r — 1) in the key identity (44) and using 
A — — j(0), we can derive 

{3{r - 1))% = (^ej: {j{s)mrLo - m)% + o{e'). (6?) 

s=0 

In the case t — 1, this expression yields 

{m)%-lPE{j{oryi^+o{E^). (68) 
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Substituting this into (67), we obtain 

(j(r - 1))% = PEJ:9{s) (j(«)j(0))Lo + 0{E'), (69) 



s=0 



where 9(s) = 1 for s > 1, and ^^(0) = 1/2. When we compare (66) and (69), 
using the time correlation function 



c(t) ^ mmrLo , m 

we find the relation for t > 0: 

e{t)C{t)^TR{t). (71) 

(71) can be rewritten as 

C{t) = T{R{t) + R{-t)) (72) 

for all t. The relation (72) is called the fluctuation-response relation, which is 
a dynamic extension of the Green-Kubo relation. 



5.2 Main claim (time- dependent potential) 



The relation between the time correlation function and the response function 
can also be derived when a perturbation potential function becomes time 
dependent. Specifically, let us consider a situation when a time-dependent 
potential 

T4 = v{t) sin -J- (73) 

is applied to the system. That is, the equilibrium Hamiltonian (16) is modified 
as 

L 

Hv{'n,t)=Ho{r^) + Y.^-y-- (74) 

x-=l 

Here, the amplitude of the density fluctuation with wavenumber 27r/L is 

^ 2tix 
m = T.Vx{Lt)sm--. (75) 

x=l ^ 

Using (75), we define the response function, Rp{t), by the relation 

oo 

m) = -i:Rpi^Ht-^)+o{A'), (76) 

s=0 

with the initial condition given at t = —oo. Note that (p(t)) denotes the 
average in the situation. Rp{t) — ior t < due to the causality. Then, for 
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t > 0, we can derive the relation 

C,it)=TR,it), (77) 

where 

The relation (77) is also called the fluctuation-response relation. 
We provide the proof of (77) using a special situation 

v{t) = A for t > 0, 

= for t<0. (79) 

In this case, the relation (77) is explicitly rewritten as 

where ( )y denotes the statistical average by the canonical distribution of 
V = at t = and by the transition probability T with w = A for t > 0. 
(Note that the canonical distribution is realized at t = 0, because the initial 
condition is given a,t t — — oo.) The form (80) can also be derived using a 
method similar to that used for the derivation of the Green-Kubo relation. 
Indeed, by repeating a procedure similar to those in Subsection 4.1, wc can 
obtain the following key identity for an arbitrary history dependent quantity, 

(A)*^ = (e^^('^«-^(o))i)*^. (81) 

Setting A([t7]) = pit) — p(0) and expanding the right hand side with respect 
to A, we obtain (80). 

Remark: Although (72) and (77) have similar forms, it should be noted that 
Rp does not represent the current response to the time-dependent potential 
perturbation. When we consider the current response to the timc-dcpcndcnt 
potential perturbation, we define the response function of the time difference 
ofp, i?d(t),by 

oo 

{pit + 1) - m) = - E Rii^Wt - + O(A^). (82) 

s=0 

Note that the time difference of p is related to the current due to the continuum 
equation of the density field. Using the response function, Rp{t), the relation 
(77) can be written as 

C,(t + l)-C,(t) = Ti?J(t). (83) 
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Remark2: In addition to the proof using the key identity (81), which is directly 
due to the local detailed balance condition, there is another method to prove 
(80) [40]. Let Py"(T7) be the canonical distribution of the system with the 
perturbation T4 of v{t) — A. It is easy to find 

pT{rj) = (1 - /3Ap) p^^^T)) + 0(A2). (84) 

Using this, we obtain 

m - m)v = am - m) (i + ^ApmrUv + o(a') 

= /3A mm)lv=o - /3A {P%y^, + OiA'). (85) 

(85) is equivalent to (80). In this method, we only use the condition that the 
stationary distribution of the system with the potential modulation is writ- 
ten as the canonical distribution. It means that (80) can be derived without 
using the local detailed balance condition, unhke the case of the Green-Kubo 
relation. 

5.3 Fluctuation relation 

When the limit t — > cxo is taken, (80) becomes 

(P>Lo,. = 4(''')to,..o + 0(A^). (86) 

As shown in the above proof {Remark2), this kind of relation generally holds 
due to equilibrium statistical mechanics. The fluctuation-response relation 
(80) is a dynamic extension of this relation. 

Here, we emphasize that (86) can be understood from a thermodynamic point 
of view. When L is sufficiently large, Vx of v{t) = A defined by (73), slowly 
varies in x. Thus, and {tIx)%^qv regarded as smooth functions V{x) 
and p{x). Then, from a thermodynamic variational principle and using the 
free energy density for the equilibrium lattice gas, /(T, p), the most probable 
density profile p{x) under the influence of the slowly varying potential V{x) 
is obtained as the minimizer of the functional 

Fm^)}) = r dx[f{T,(t>{x)) + V{x)(t>{x)l (87) 

under the constraint condition 

dxMx) = Lp. (88) 

Jo 
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Prom this variational principle, p{x) is determined so as to satisfy 

^i{T,p{x)) + V{x)-X = 0, (89) 

where A is the Lagrange multiphcr due to the constraint condition, and ii{T, p) 
is the chemical potential defined by 

Difi^erentiating the relation (89) with respect to x, we obtain 



dp{x) dp{T,p) 



, ^ dX 

p=p{x) 



dx dp 

Substituting (73) into (91), and using (75), we derive 



+ = (91) 



(P>S,. = -^(^1 +0(A^). (92, 



where p is the averaged density. 

Now, let us recall that the thermodynamic variational principle is derived from 
the Einstein-Boltzmann formula for the probability of the density profile: 

Prob({p(a;)}) ~ 5 dx{p{x) - p)j e-''^«''(=^)» . (93) 

When V{x) — 0, setting p{x) = p + Sp{x) and expanding f{T,p{x)) with 
respect to Sp{x), we obtain 

Prob({5p(x)}) ~ S (^£dx5p(x)j e-i(^^)/o^'^^^'^^^^^)', (94) 

using the definition of the chemical potential (90). The expression (94) leads 
to the fiuctuation relation 

where x is the intensity of density fluctuations deflned by 



st 

X^Ax{(5p(x)r)^^^ (96) 

in the continuum description. Here, Ax is the length much longer than the 
microscopic correlation length, but much shorter than the system size. p{x) 
represents the coarse-grained density over the region of the length Ax. 
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Substituting this result into (92), we obtain 

(Pri=o,y = -^X + 0(A^). (97) 
Thus, (86) is the fluctuation relation, if we can prove the relation 



Proof of (98): Let us define the intensity of the density fiuctuations x more 
explicitly than that in the continuum description given in (96). Consider a 
coarse-grained density of the subsystem with the size t. 

Pi = )Y.^.- (99) 

^ x=l 

Then, using the quantity 

d,,^ = £((p,-p)2)'^^^, (100) 
the intensity of density fluctuations, x, should be deflned as 

X = lim lim di^i- (101) 



Based on this deflnition, we now derive (98). Notice the asymptotic relation 

di,L X-^ (102) 

for 1 <^ £ <^ L. The derivation of (102) is explained as follows. 

Let Ni = pii and Ni = pi. Then, according to the central limit theorem, the 
probability of approaches to the form 



(N^-N^)"^ (N-Nf-jN-Nf))"^ 

Prob(A^^) ~ e 2x(L-e) 

= e ^ (103) 

when i and L — i become large. Prom this, the probability of pe is written as 

Prob(p^) ~ e ^x(L-e) . (104) 

This implies (102). 
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Using this asymptotic relation, we obtain 



X = 2 lim dL/2,L- 



(105) 



Furthermore, we can derive 




St 



E=0 ' 



(106) 



using a straightforward calculation. Combining (105) and (106), we obtain the 
expression (98). 



6 Einstein relation 

As is well-known, one of the celebrated papers written by Einstein in 1905 
is regarded as genesis of linear response theory [41,42]. Even in many-body 
systems, if a diffusion constant can be properly defined, it is connected with the 
conductivity in an extended form of the original Einstein relation for one-body 
systems. In this section, we derive the Einstein relation for the one-dimensional 
driven lattice gas. We also explain the interrelations among the Grccn-Kubo 
relation, the fluctuation-response relations, and the Einstein relation. At the 
end of this section, we discuss the Einstein relation from the thermodynamic 
point of view. 

6.1 Main claim 

We begin by defining the density diffusion constant, D. We first assume that 
the most probable process of the density field is described by a diffusion equa- 
tion. On this assumption, D should be defined as a coefficient of the diffusion 
term that appears in the equation. As an example, we consider the time evo- 
lution of the density field after the slowly varying weak potential is turned on 
at t = 0. When the most probable behavior in this situation can be described 
by the diffusion equation with the diffusion constant D, the response function 
Rp{t), given in (76), can be expressed as 



In this paper, we assume that (107) determines the value of D. 

Now, the Einstein relation in this model connects the diffusion constant, D, 
with the conductivity, a, which is defined by (37), as 




(107) 



Dx = aT. 



(108) 
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We derive (108) using both the Green-Kubo relation (40) and the fluctuation- 
response relation (80). 



Proof: From (80) and (107), we obtain 

((p(*)-p(0))')-;^„ = 2(p^);^Jl-e-"tt)'']. (109) 

Notice the continuity equation 

r]^{Lt) - 77^(0) = - J^,t + J^_i,t, (110) 

with 

J.,t = 'Z^.{ri{k-l)^7j{k)), (111) 

k=l 

where ^xi^l v') is given in (28). Using (110) and (111), we derive 

m - m = (^) E J^,t cos ^ + o (-1) . (112) 

Next, we define 

= lim lim Bi,t(L) (113) 

t— »oo L— >oo 

with 

where the order of the two limits in (113) must not be exchanged. Then, from 
(98), (109) and (112), we find 



2 / o \ St 

= Dx. (115) 

Furthermore, we can derive B^o — B because the {Jx,tJy,t)%=o is a rapidly 
decaying function of \x — y\ (<^ L/2) at a fixed t. Thus, we obtain 

B = Dx. (116) 

Combining this relation with the Green-Kubo relation (40), we arrive at the 
Einstein relation (108). 

Remark: We comment on the diffusion of a marked particle which is character- 
ized by a tracer diffusion constant. The tracer diffusion constant is not equal 
to the bulk diffusion constant, D, represented by (107), and it is not related 
to the other quantities presented in this paper. 
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6.2 Interrelation 

We argue the interrelations among the Green-Kubo relation (40), the fluctuation- 
response relation (80), and the Einstein relation (108). To simphfy the argu- 
ment, we assume that the time correlation function, Cp{t), is written in an 
exponential form as a function of time. That is. 



Based on this assumption, the fluctuation-response relation (80) in the simple 
form is 



from the continuity equation. Now, when the fluctuation-response relation 
in the form (118) is vahd, we find with the aid of (119) that the Green- 
Kubo relation (40) and the Einstein relation (108) are equivalent. Also, when 
both the relations (40) and (108) are valid, we obtain the fiuctuation-response 
relation (118) by using (119). In this way, we have confirmed that any two of 
the three relations lead to the other relation, when Rp{t) and Cp{t) have the 
exponentially decaying forms. 

For the driven lattice gas, using the diffusion constant, D^, given in (117), 
Katz et al presented the rigorous proofs of the Green-Kubo relation and the 
Einstein relation [21] (although in a less systematic way than that presented 
here). It should be noted that there is the fluctuation-response relation (118) 
behind the adoption of this definition of the diffusion constant. Related to 
this issue, it is problematic in a system far from equilibrium whether to define 
the diffusion constant from a response function or a time correlation function, 
because for the one-dimensional driven lattice gas, the fluctuation-response 
relation (77) is violated in the nonequilibrium steady states [9]. 

Below, we demonstrate another important interrelation; the Einstein relation 
is equivalent to the thermodynamic fluctuation relation (95), when the re- 
laxation of the density field is described more phenomenologically instead of 
using (107). 

Let J(p, E) be the averaged particle current in a system with density p and 
driving force E. When the slowly varying potential \4 of v(t) = A, given 
in (73), is applied to the system of ii^ = 0, first the current should become 
J{{'nxYE=Qy , -dV^/dx). Consequently, a diffusive current D {{rj^+i - r/a))^=o,y 




(117) 



Similar to the derivation of (115), we can derive 



(118) 



(119) 
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should appear so that the total current is equal to zero. That is, 



J{{Vx)t=0,V ' - ^ ((^7^+1 - Vx)TE=oy = 0. (120) 

Precisely speaking, this relation is assumed to define D phenomenologically. 
As explained in Section 5, {rix)%=Q y and T4 can be regarded as p{x) and V{x). 
Thus, (120) can be rewritten as 

The expansion of (121) with respect to dV/dx yields 

-a^-iP^ = 0. (122) 
dx dx 

Note that \dV/dx\ is sufficiently small because V{x) slowly varies in x. Substi- 
tuting the result of the thermodynamic variational principle (89) into (122), 
we obtain 

a^-^-D^^. (123) 

op 

Then, the thermodynamic fiuctuation relation (95) leads to the Einstein rela- 
tion (108). 

The above argument states that the non-triviality of the Einstein relation 
depends on the selection of which D to use among (107), (117) and (121). 
Also, it is amazing to sec that the thermodynamic variational principle, the 
thermodynamic fiuctuation relation, and linear response theory are mutually 
related. 

Here, let us recall that the fiuctuation-response relation (77) can be derived 
without the local detailed balance condition, but only if a general Markov 
process compatible with the statistical mechanics is used. Using this fact, and 
admitting the phenomenological derivation of the Einstein relation, we can 
derive the Green-Kubo relation. It implies that all the relations in linear re- 
sponse theory can be understood without the local detailed balance condition. 
Although the argument is correct, it assumes that the density field obeys a 
linear diffusion equation. We wish to emphasize that the Green-Kubo rela- 
tion and the fluctuation-response relation are derived using the local detailed 
balance condition even when the density field does not obey a simple linear 
equation. 
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7 Green-Kubo relation in the boundeiry-driven lattice gas 



In this section, we study a different type of nonequihbrium lattice gases where 
the particles are not driven by an external driving force in the bulk, but are 
driven by the difference in chemical potential between two particle reservoirs 
in contact with the boundaries. Due to the local detailed balance condition 
imposed on the model, we can derive the Green-Kubo relation in a boundary- 
driven lattice gas. 

7.1 One- dimensional boundary- driven lattice gas 

Consider the system in which the sites x = 1 and x = L are in contact with the 
particle reservoirs that have the chemical potentials /xi and /^l- Specifically, 
we set III = II and a*l = A* + A/x. When A/x ^ 0, the particles flow in one 
direction on average, but when A// = 0, the grand-canonical distribution is 
realized. 

We describe the time evolution of the configuration of the system, r], by iter- 
ating the following three Monte Carlo steps: 

In the first step, a particle is created and annihilated at the site x = 1. This 
process is represented by hipping the value of rji with the probability Ci{r}). 
Here, changing from 771 = to 771 = 1 corresponds to creating a particle, 
while changing from 771 = 1 to 771 = corresponds to annihilating a particle. 
The flipping probability, ci{r}), must satisfy the condition, which is a simple 
extension of (17), 



where Ho{r]) represents the interaction between particles, and its form is given 
by (16). Recall that ij^ denotes the conflguration of the system after flipping 
the value of 771 (see Section 3 for the usage ). Logically speaking, the condition 
(124) is necessary for the transition probability to satisfy the local detailed 
balance condition in nonequihbrium steady states. We will check this condition 
in the next subsection. 

In the second step, we provide a time evolution rule for the bulk. We randomly 
select a nearest neighbor pair x, y, and exchange the value of 77^; and 77^; with 
the probability Cbuik(2^, V, v) to satisfy 



We also deflne Cbuik(a;, yiv) =0 fo^ \^ ~ y\ 1- Note that this process is the 
same as that for the lattice gas without the particle reservoirs. 




(^1 ) g-/3[J^o (rji ) -i?o (rj)- Ml (C'?' )i -m )] ^ 



(124) 



Cbuik(x, y; v) = Cbu,k(x, y; rj^^) e-/^[^o(^-)-/^o . 



(125) 
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In the third step, the procedure of the first step is performed at the site x = L. 
That is, we change the value of r]L to that of 1 — 77/, with the probabihty cl(»7) 
to satisfy 

CLirj) = CL(77^)e-''[^°(^'')-f^o(^)-'*i'«^^)^-''^)]. (126) 

Then, the corresponding transition probabihties Ti{ri ry'), Tbuik(^ — ^ ij') 
and, Tl(t7 — > 77') are defined as 



T.irj rj') = c^ir1) for rj' = rj' ^ rj, (127) 
7buik(r7^r7') = -^Cbuik(x,y;r7) for -n'^r^-y^n, (128) 
7l(r7 ^ V) = CL(r7) for 77' = 77^ ^ r/. (129) 

The transition probabihties Ti(77 77), Tbuik(^ 11)1 ^"^^ Ti^ijl ~^ v) 
determined from the normahzation conditions of the probabihty 



5:ri(77^r7') = l, (130) 

r,' 

E^buik(r7^r7') = l, (131) 
n' 

E^l(^-V) = 1, (132) 



respectively. Notice that these expressions of the transition probabilities lead 
to the equalities 



Jbuik(r? -» r?0 ^ piHo{v')-Ho{v)] ^34) 



7.2 Local detailed balance condition 



Using the above time evolution rule, a history of configurations during r MCS 
is denoted as [77] = (77(0), 77(1), ■■ • ,77(3Lr)). For a segment of the history 
[v]k = (rtiSk), 77(3/c + 1), 77(3/? + 2), 77(3/^ + 3)), where = 0, 1, 2, • • • , Lr - 1, 
we define the product of the transition probability, T([77]fc), as 
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T{[ri]k) = TiirjiSk) ^ Ti{3k + l))Tbuik(r7(3A; + 1) ^ rjiSk + 2)) (136) 
Associated with T{[r]]k), we also define 



T*i[v]k) = TL{r]{3k) ^ r]{3k + l))Tbuik(r7(3fc + 1) ^ ri{3k + 2)) (137) 
xTi{'n{3k + 2) ^rj{3k + 3)). 

The transition probabihty T*{[r}]k) corresponds to the time evolution rule ob- 
tained by exchanging the first and third steps introduced in the last subsection. 
Furthermore, [77]^ = (r7(3A; + 3), t/(3A; + 2), r7(3/i; + l), t/(3/c)) represents the time 
reversed segment of [77]^. Then, from (133), (134) and (135), we obtain the 
relation 

^ g-/3[i/o(r,(3fc+3))-/fo(»7(3fc))-Mi*i(fc)-ML*L(fc)] QogN 

T*{[il]l) ' ^ ^ 

where 

^,{k)=r]i{3k + l)-ri,{3k), (139) 

and 

q>^{k)=riLi3k + 3)-riLi3k + 2). (140) 

$i(A;) and ^h{k) are numbers of particles that fiow into the system from the 
particle reservoirs during the time steps from 3/c to 3/c + 3 at a; = 1 and at 
X = L, respectively. 

Now, the increment of the particle number at the site x from t — 3k + 3 to 
t = 3k is 

5ri,{k)=r,,{3k + 3)-r,^(3k). (141) 
Because the total particle number is conserved, the following relation holds: 

L 

Y,5Vx{k) = ^i{k) + ^i^{k). (142) 

x=l 

Using this relation, and recalling that /ii — /i and //l = A* + A/x, the equality 
(138) can be rewritten as 



-(3[Ho {n{3k+3))~Ho {vm)-fM (fe)- A^^l (k)] 



^e-^[(^o(''(3fc+3))-^^f^j7,,(3fe+3))-(/fo{rj(3fc))-^^f^jr,,(3fe))-A^*L(fe)] _ ^^^3^ 

This is regarded as the local detailed balance condition in the boundary- driven 
lattice gas. 

In the equihbrium case A// = 0, (143) leads to 

p^--(^(3A;))T([r7]fc) = p^^-'^{ri{3k + 3))T*([r7]*), (144) 
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with the grand-canonical distribution p^^^^{r]) 

^gcan^^^ = - e-z^I^oW-z^ELi"-], (145) 

z 

where Z denotes the normahzation factor. The equahty (144) corresponds to 
the detailed balance condition with respect to the grand-canonical distribu- 
tion. This ensures that (145) is the steady state distribution of the stochastic 
model. 



7.3 Main claim 



In order to derive the Green-Kubo relation in this model, we assume that the 
grand-canonical distribution ^§'=^'1(77) of the chemical potential /i is realized at 
t — 0. Then, we consider the situation where A/x 7^ when t >0.ln this case, 
the probability of a history [rj] — (t7(0), rj{l), • • • , rj{3LT)) is given by 

PaM) ^P^''"\rimn[ri]o)---T{[riUr-i). (146) 
The path probability associated with T* is also defined as 

PrM) =P'^"Xr7(0))T*([r7]o) • ■■T*{[nU_^). (147) 
Noting the equality 

^a;([^]) = P^^^'Xr7(^r))T*([r7]L-i) • ■■T*{[riU (148) 
and using the local detail balance condition (143), we obtain 



with 

Miv]) = - E ^l(^). (150) 



Lt-1 



k=0 



Then, we write the statistical average in the transient process by the proba- 
bihty, P'^M)^ ^ 

{A)t^j:PAM)M[ri]), (151) 

iv] 

for a history dependent quantity ^([ry]). (^)^^ is also defined in the same way, 
replacing P^i,{[rj]) by Pam(W) in (1^1) • 

Noting the equality Jr ([*?]) = "-'r ([*?]), we obtain the key identity 
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Iv] 
M 

^-pA^rJAM) a)''*, (152) 



where we have defined A([ry]) = ^([77]). 

Here, setting A = J^- in (152), taking the hmit r — > 00 and expanding the 
right hand side of (152) with respect to A/i, we can derive 

{Jr)t, = ^^H + 0{An'), (153) 

with 



B = hm -{(Jr) ) ■ (154) 

Note that ( and ( )^* arc equal to the average in the steady state ( 
with the large r limit and a fixed system size L. The expression (153) pro- 
vides the Green-Kubo relation in the boundary- driven lattice gas, because the 
current correlation B in the system under the equilibrium condition is related 
to the conductivity a in the nonequilibrium in the form B — aT, where a is 
defined as 

<7^1imi4^. (155) 

We remark that the other relations derived in Sections 5 and 6 can be obtained, 
using the local detailed balance condition (143). As far as we know, the Green- 
Kubo relation in this model was presented in a different manner in Ref. [20]. 
Our derivation method is more pedagogical in the sense that the role of the 
local detailed balance condition is explicit. 



8 Discussion 



We have derived the universal relations of linear response theory for bulk- 
driven and boundary-driven lattice gases in a simple manner. We have also 
elucidated the interrelations among these universal relations, statistical me- 
chanics and thermodynamics. However, in the arguments presented above, 
some important topics related to linear response theory were not discussed. 
Among them, there are two that are particularly worth consideration, the 
derivation of the universal relations from classical or quantum systems, and 
the extension of the universal relations to forms valid in nonequilibrium steady 
states far from equilibrium. In this final section, we present remarks on these 
two topics. 
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8.1 Microscopic understanding 



It is well known that the universal relations of linear response theory can be 
derived from a microscopic description in the following way. (i) The system 
under consideration is assumed to be in equilibrium at t — — oo; that is, the 
distribTition of microscopic variables is assumed to be canonical at t = — oo. (ii) 
An external force is applied to the system under the assumption that the time 
evolution obeys the Liouville and von Neumann equations. This assumption 
implies that the system does not interact with other dynamical degrees of 
freedom, (iii) The statistical average of the quantity of interest at time t is 
calculated perturbativelly up to linear order in the applied external force. 

Within such a framework, it is straightforward to derive the hnear response 
relations, but it is difficult to understand the physical picture. For example, 
let us consider an electric conduction system. In such a system, energy is 
continuously introduced into the system when an electric field is applied. Thus, 
obviously, a steady state is never realized unless the system contact with a heat 
bath. Then, recalling the step (ii) described above in the previous paragraph, 
it can be claimed that the formal calculation is not performed based on the 
realization of nonequilibrium steady states, even if the desired formulas can 
be derived by this method. 

Furthermore, it should be noted that the time correlation (e.g. current corre- 
lations) in the formal perturbation scheme is calculated by use of mechanical 
time evolution equations without identifying the degrees of freedom that con- 
stitute the heat bath. If the energy dissipation into the heat bath is taken into 
account when the time correlation is calculated, it is plausible that the result 
of the calculation would depend on the nature of the heat bath. Indeed, in 
the stochastic models studied in this paper, the functional form of the time 
correlation depends on the choice of the stochastic rule (e.g. (19), (20), and 
(21)), which represent the nature of the heat bath, though the validity of the 
universal relations is independent of this choice. 

Considering these points, we find that we do not obtain a clear understand- 
ing of nonequilibrium systems, even in the linear response regime, when we 
study systems on the basis of a microscopic description. Here, let us recall 
our method of derivation for the universal relations in the stochastic models. 
In our derivations, the steady state is prepared from the outset, and all the 
relations can be obtained by use of the local detailed balance condition. Thus, 
focusing on the realization of nonequilibrium steady states and the local de- 
tailed balance condition, we suggest that the two problems described below 
should be studied in microscopic systems. 

The first problem is to clarify the conditions under which a subsystem exhibit- 
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ing a noncquilibrium steady state is determined in a purely mechanical system 
despite a continuous injection of energy. (Note that a nonequihbrium system 
with a Gaussian thermostat [43] is not regarded as a purely mechanical system 
of the type we wish to study) Then, we characterize the heat bath by studying 
the mechanism of the energy dissipation, but it is not certain that a simple 
characterization is possible. If wc succeed in the proper characterization, the 
second problem is to demonstrate the vahdity of the local detailed balance 
condition for these states. Then, given the local detailed balance condition, 
it is straightforward to derive universal relations. In this way, through the 
present study, we have identified a new direction in the investigation of the 
nature of nonequihbrium systems on the basis of a microscopic description. 



8.2 Extension of the linear response relations 

One may wonder whether it is possible to extend linear response theory to 
systems far from equilibrium. A partial answer to this question is provided by 
the nonlinear response relation given in (44), which provides a starting point to 
connect transport properties of systems far from equilibrium with dynamical 
fiuctuations. Indeed, when Kawasaki and Gunton derived an expression for the 
nonlinear shear viscosity from the Liouville equation, the nonlinear response 
formula [26], which takes essentially the same form as (44), played a key role. 

There is another approach to the extension of the universal relations that does 
not involve nonlinear response theory. In this approach, instead of studying 
nonlinear transport properties, the linear response properties near nonequi- 
hbrium steady states are studied by applying small perturbations as probes. 
It is believed that there is a close relationship between the linear response 
properties and the dynamical fiuctuation properties, although, obviously, the 
universal relations that hold within the linear response regime cannot be valid 
in general outside this regime. Unlike the nonlinear response theory repre- 
sented by (44), this type of extension of the universal relations beyond the 
linear response regime has not yet been established. 

As a preliminary step in the attempt to employ this second type of approach 
to construct an extension of the universal relations to systems far from equilib- 
rium (i.e. outside the linear response regime), recently, numerical experiments 
on a two-dimensional driven lattice gas were performed [11]. According to the 
results of this study, with respect to the properties of the system along the 
direction transverse to the external driving force, the linear responses of the 
system to perturbation forces are related to the fluctuations in the same form 
as the linear response relations. Here, it is important to note that in the direc- 
tion transverse to the external driving force, the thermodynamic free energy, 
from which we can obtain the probability distribution for density fiuctuations. 
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was constructed using the Maxwell relation [44,45]. (The Maxwell relation is 
the integrability condition yielding the free energy.) We conjecture that the 
validity of the extended forms of the linear response relations [11] is related to 
the existence of the thermodynamic function constructed in Ref . [44] , recaUing 
the argument given in Section 6. 

Although these numerical results for states far from equilibrium are interest- 
ing, our understanding of the extended relations remains poor. In particular, 
there are the following two serious problems. First, it is well known that long- 
range spatial correlations of density fluctuations exist in nonequilibrium sys- 
tems of two or more dimensional [46,47,48]. These long-range correlations are 
inconsistent with the extensive nature of thermodynamic fluctuation theory. 
In fact, we do not understand how the numerical results of Ref. [11,44] can 
be reconciled with the existence of long-range correlations. (See Ref. [45] for 
a related discussion). Second, while we were able to numerically construct an 
extended free energy and conflrm the validity of the extended linear response 
relations along the direction transverse to the external driving force, it seems 
more difficult to find similar relations along the direction parallel to the driv- 
ing force. Indeed, the study of a one-dimensional driven lattice gas far from 
equilibrium has revealed that the relationship between the fluctuations and 
the response to a probe force is complicated, because the fluctuations take the 
influence of hydrodynamic effects [9] . 

In order to construct a self-contained theory of nonequilibrium steady states, 
we need a unifled treatment of thermodynamic and hydrodynamic fluctua- 
tions. For example, thinking optimistically, despite the existence of the long- 
range spatial correlations in nonequilibrium steady states, it might be possible 
to extract the thermodynamic component of the fluctuations by separating out 
of the long-range component. At present, however, there is no theory providing 
such a treatment. However, we should mention that the additivity principle 
proposed in Refs. [49,50] might provide a useful tool to study this problem, be- 
cause this principle provides an elegant method by which the thermodynamic 
part of density fluctuations can be distinguished from the hydrodynamic part. 
Although the validity of this principle has been confirmed only for certain 
exactly solvable nonequilibrium models, it would be interesting to consider its 
application to a wider class of nonequilibrium models and to study how ther- 
modynamic fluctuations can be extracted by use of this additivity principle. 

If we could construct a unified treatment of fiuctuations for nonequilibrium 
steady sates far from equilibrium, our next step would be to seek new univer- 
sal relations between fluctuations and response properties to probe forces in 
such systems. In that pursuit, and in the construction of a general theory of 
nonequilibrium steady states, the present paper, in which we have elucidated 
the interrelations among the linear response relations, statistical mechanics 
and thermodynamics, should serve as a useful guide. 
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